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ABSTRACT 


This dissertation introduces the idea of equicontrollability and 
studies its application to the linear time-invariant model-following 
problem. The problem is presented in the form of two systems: generically 

called the plant and the model. The requirement is to find a controller 
to apply to the plant so that the resultant compensated plant behaves, 
in an input- output sense, the same as the model. All systems are 
assumed to be linear and time-invariant. 

The basic approach used is to find suitable equicontrollable 
realizations of the plant and model and to utilize feedback so as to 
produce a controller of minimal state dimension. The concept of 
equicontrollability (introduced here) is a generalization of control 
canonical (phase variable) form applied to multivariable systems. 

It allows one to visualize clearly the effects of feedback and to 
pinpoint the parameters of a multivariable system which are invariant 
under feedback. 

The basic contributions contained in this work are: (1) 
the development of equicontrollable form; (2) solution of the model- 
following problem in an entirely algorithmic way, suitable for computer 
programming, and (3) resolution of some questions on system decoupling 
(along with the application of the above algorithm to accomplish 
decoupling, as shown in Appendix C). 
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I . INTRODUCTION 


A. Description of the Problem 

The model-following problem has appeared periodically in the 
literature for some time, it most often is referenced in conjunction 
with proposed solutions for flight control problems. These have the 
property of having no a priori cost function to be minimized (contrary 
to most problems in modern control theory). Instead, a model is 
specified whose dynamic response is considered desirable. The problem 
then becomes finding a controller (compensator) which, when added to 
a given plant, will cause the resultant system to have a response as 
close to the model's as possible. 

A typical example of this kind of problem might arise during 
the design of an SST. The future pilots of the SST would like to fly 
the aircraft before it is built. This clear contradiction is generally 
"solved" by construction of a simulator which does a credible job of 
reproducing the "feel" of the aircraft, but is still not quite the real 
thing. A recent proposal has been to take a small jet transport and 
build a complete SST nose section on its front. Then some sort of 
artificial feel system would connect the controls to the transport in 
such a way as to make it feel to the pilots as if it were an SST 
behind them. In this problem the model is the SST, the plant is the 
transport and the artificial feel system is the compensator to be 
designed. This problem was considered by Rynaski and Whitbeck [17] 
with some success. 

Another problem of the same type occurs in the design of 
compensation for VTOL aircraft. In this case the plant is the helicopter 
or VTOL and the model is a mathematical description of what sort of 
system the pilot would like to fly. Such a description is probably no 
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more detailed than the requirements of asymptotic stability, negligible 
overshoot and decoupling of lateral and longitudinal dynamics. Time 
constants are stated equally as grossly. Nonetheless, a model may be 
formulated which the pilot would accept. This type of problem has been 
considered by Wolovich and Shirley [21]. Unfortunately, they had to 
compromise the model substantially to get a solution. 

In summary then, the model-following problem involves some sort 
of plant and a model whose response the compensated plant is to emulate. 
The key factor in every problem statement of this type is a lack of a 
measure of performance. No cost function is given to minimize, but 
rather the designer wants the model and compensated plant to be the 
same in response; or, barring that, "as close to the same as possible". 
This is the point on which the problem hinges. At such time as the 
measure of error is defined, the problem is half solved. 

One might notice that there is a definite link between the 
model-following problem and classical control theory; at least in so 
far as the observation that virtually all classical design problems 
specify enough desired parameters to roughly define a model. For 
example, such criteria as overshoot, risetime, gain margin and so 
forth can be translated into the specification of a second or third 
order system. Certainly no cost function is given to be minimized 
and the designer would like the plant and such a model to behave as 
nearly the same as possible. Hence, the original form of the model- 
following problem acts as a sort of bridge between what is commonly 
called optimal control, i.e. problems with a cost functional, and 
classical control problems. 

Strangely enough, this link was to some degree born out 
historically. At the time when a type of model-following problem was 
being considered in a frequency domain sense (see Freeman [8], 

Kavanaugh [11] and Morgan [ 14]) ."modern" control theory was waxing; 
in particular Kalman’s work was making the quadratic performance index 
de riguer . The net result was a loss of the link between the two 
approaches. Later some authors attempted to re-establish it (notably 
Brockett [3]), but most research continued in the direction of optimal 
control theory. 
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B. 


Proposed Solutions 


1. The Frequency Doma in Approach 

The initial interest in the model-following problem came from 
the relative ease with which its scalar version could be solved. That 
version involved the block diagram in Fig. 1.1: 


r(t) 



y(t) 


Fig. 1.1 

The Scalar Problem 


r(t) is a scalar reference input, V(s) is a given linear model and 
W(s) is the plant. The plant has been imbedded into a unity feedback 
loop with the series compensator C(s) . The objective is to find 
C(s) so that y(t) = p(t) . In other words so that: 


C(s)W(s) 

1 + C(s)W(s) 


= V(s) 


or that: 


C(s) = 


_1 # V(s) 

W(s) ’ l-V(s) 


( 1 . 1 ) 


There were, of course, some restrictions placed on C(s) . Most often 
they would be: 

1. the closed loop system using C(s) must be stable and, 

2. C(s) must be realizable (or, more restrictively , low pass). 

These two restrictions are easily dispatched. First, if the numerator 
of W(s) has a root in the RHP, then clearly we must require that V(s) 
also have one at the same location. This assumption will guarantee 
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stability. The second restriction is no more difficult to deal with. 

If the numerator of C(s) is of greater dimension than its denominator, 
one may modify the model in the form a r V(s)/(s+a) r yielding: 


C(s) - — 
W(s) 


a r v(s) 

(s+a) r -a r v(s) 


( 1 . 2 ) 


We may certainly choose an r so that C(s) will be realizable. Then 
a can be chosen so that the effect on the model will be as small as 
desired (at the cost of increased feedback gain of course). 

The manipulations engendered by the scalar case are fairly 
transparent, while the multivariable case is quite opaque. The block 
diagram may be left the same, except that all variables would be 
vectors instead of scalars, and all transfer functions would be matrices. 

We may then formally derive the required compensator: 

C(s) = [W(s)] -1 V(s)[I-V(s)] _1 (1.3) 

assuming that all inverses exist. Now, unfortunately, we cannot so 
glibly arrange to satisfy the two given requirements. Stability is not 
so clear, and realizability is more of a chore. Worse yet, the computation 
of the inverses is especially complicated. All in all, the problem, as 
given, is very difficult. It was only saved from extinction by the 
advent of the quadratic loss approach which circumvented the difficulties 
via a reformulation. Although we shall also reformulate it, Chapter V 
contains a consideration of a similar problem which the algorithm given 
in the sequel can solve. 

2. The Time Domain Approach 

The bulk of work in attempting to solve the model following 
problem has used optimal control theory. To do this requires the 
generation of a cost function. Since the actual problem statement does 
not include one, it is necessary for the designer to create one to suit 
both his convenience and to produce an acceptable answer. Most all of 
the papers written in the area specify a quadratic performance measure. 

More specifically, if the stateof the model is | and that of the plant 
is x , then thsy chcossi 
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J = J ^(x-£) T Q(x-!) + U T RU j dt (1.4) 

where u is the control applied, as their measure. Since the plant and 
model are assumed to be linear constant systems, the solution to such a 
problem is well known. Such an approach has been used by Kalman and 
Englar [10], Rynaski and Whitbeck [16], Rynaski at al_ [17], Tyler [ 18] , 
Tyler and Tuteur [19], Asseo [1], and Winsor and Roy [20]. These have 
been, in some sense, ordered as to increasing sophistication in specify- 
ing the matrices Q and R ; i.e.: the relative cost or states and 

control effort. Coincidentally, they also appear in nearly correct 
chronological order and order of increasingly useful results .(a judgment 
based on the author's personal bias and minimal data). 

The key difficulty in all of these attempts is the lack of 
definition of Q and R . This is combatted by methods ranging from 
cut and try to the addition of side constraints. None is especially 
preferable to another, except that perhaps Winsor and Roy's approach 
appears to be slightly more of a science than an art. Therein lies the 
biggest fault in any of them. The production of the control law is 
more of an art depending on the wit, cunning and experience of the 
engineer rather than an algorithm which a computer can be taught to solve. 

At the same time, the very nature of the solution lends itself 
to criticism. Those familiar with the solution to the linear quadratic 
loss regulator problem recall that it consists of constant feedback from 
all the states. Sipce the designer here had to use the states of both 
the system and the model in the cost function, he must build a realization 
of the model as part of the controller. Diagramatically , the compensated 
plant is shown in Fig. 1.2: 
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Fig. 1.2 

Compensated Plant 

The K would be produced by the solution to the quadratic loss problem 
created. Note that will not be zero, even in the case where a 

K 2 exists such that the system in Fig. 1.3 



Fig. 1.3 

Compensated Plant, Trivial Case 

has the same response as the model I Similarly, intermediate solutions, 
such as a compensator needing only one state of dynamics to do the job 
perfectly, will not be uncovered. In summary, the solutions of this 
type suffer these difficulties: 

1. Since the cost function is unspecified, it must be 

chosen ad hoc . Moreover, since the complete solution to 
the inverse problem is as yet unknown, iteration toward a 
solution is often a blind search. 
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. Exact solutions, even when they exist, will not be found. 

The compensator must contain a realization of the model — 
a real annoyance. 

3. Although not mentioned here as yet, the solutions achieved 
by the above authors were none too good, unless the problems 
posed were virtually trivial. More on this below. 

Quite a radical departure from the quadratic loss approach was 
suggested by Erzberger [6], The essence of his approach is contained 
in the following observation. Of all conceivable happenstance, the 
most fortunate, from the point of view at hand, would be for a feedback 
to exist such that the plant and model would be identical. Suppose that 
the plant is governed by the vector differential equations: 

x = Fx + Gu ; y = Hx 

and the model by: 

£ = L§ + Gr ; p = H| 

Suppose there exists a K such that: 

F + GK = L (1.7) 

then we may obtain the compensated system 

x = (F + GK)x + Gr ; y = Hx (1.8) 

by letting u = Kx + r . Now: 

x = Lx + Gr ; y = Hx (1.9) 

has the same transfer function as the model. Such a K exists if 
and only if (I - GG^)(L - F) = 9*. The solution is then K = G^(L - F) 
where G^ is the penrose pseudo- inverse of G . (See [10] or [12]). 
Close observation shows that for such a K to exist, great structural 
similarity must exist between F and L . Moreover, G must have 
sufficient rank. There have been at least two such examples given in 
the literature which were solved using the quadratic cost technique [17], 

$ 

9 represents the null matrix. 


(1.5) 

( 1 . 6 ) 
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[18]. Erzberger has shown that, for this limited class of problem, 
his direct approach involves much less computation than the quadratic 
cost approach. We have extended Erzberger' s approach for use in cases 
in which the aforementioned structural correspondence is satisfied. 

Since this wcrk is not directly related to our major exposition, it is 
treated in Appendix A. 

The work of Asseo [1] is worthy of comment. Given the plant 

and model 

PLANT: x = Fx + Gu y = lx (1.10a) 

MODEL: £ = ©S + Tto p = I to (1.10b) 

the method given in [1] finds a feedback control: u = K(o + Kx + K£ , 

J v p m 

such that K and K satisfy: 
v m 

GK = T (1.11a) 

v 

G(K -K ) = $ - F (1.11b) 

m p 


where K is arbitrary. If (1.11) are satisfied, we may substitute 
P 

into (1.10a) to obtain: 


hence: 


x = Fx + G(K u> + K x + K £) (1.12) 

v p m 


(x-£) = (F + GK p )(x - |) 


(1.13) 


and K may be chosen to make the error system as stable as is desired. 
P 

(Under the usual controllability assumptions). The difficulty lies in 
trying to find a solution to (1.11). Such a solution is vital to the 
method. In [1] a canonical form (originally derived in [2]) is used as 
a tool to solve (1.11). The claim was that for any pair (F,G) , there 
exists a transformation T such that: 



re. 

n 

n-m 

" e ' 

T _1 FT = 

1 

... L . . 
1 

and T G = 

— 


- A ' 

Bj 

m 

. i _ 


(1.14) 
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In fact, such a transformation does not, in general, exist even if we 
allow: 


n-m 

, det D 4 0 (1.15) 

m 

To check that this obviates the general solution of (1.11), we need 

only observe that first defining new states as z = T x and q = T E 

p m 

solving (1.11), in the sense of [1], is equivalent to finding P 

(nonsingular), K and K„ such that: 

J. £ 

GK 1 = P _1 r (1.16a) 

GK 2 = P -1 $ - F (1.16b) 


T _1 G = 


0 


D 


Suppose 


(F,G) 


is in the form above. Then we need to solve; 


GK 


1 




r 


(1.17a) 


' e " 

= P _1 OP - 

! 

CD 

1 


r 

- K 2- 


r 

> 

td 


(1.17b) 


But to solve (1.17) requires the existence of P as in (1.14), which is 
not, in general, possible [5]. Hence (1.11) does not have a solution for 
every plant and model. This implies, unfortunately, that Theorem 2, 
part b, [1] is false. (For reference, Theorem 2, part a, is concerned 
with the case in which G is nonsingular). 

There is a further difficulty with the approach of [1]. The 
formulation given there attempts to match states. In a reasonable model- 
following problem, we would expect no net output transformation to be 
allowed. That is, if we desire to keep the error e = y-p as small as 
possible (in some sense), to consider e = y-Tp is just short of absurd. 

It is tantamount to matching, say, the pitch response of the model to the 
roll response of the aircraft. Since x = y and % = p in the formulation 
of [1], we must conclude that the plant and model transformation are the 
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same 


( T = T \ . Hence P = T = I in (1.17) and the same 

transformation must take both the plant and the model to the same 
canonical form. This further reduces the number of cases which can be 
reasonably handled by this direct approach. Overall, the approach in 
[1]» is an improvement on Erzberger's [6], but is not as general as 
needed for multivariable control. 

Winsor and Roy's paper [20] errs in description but appears 
correct in mathematics. They falsely claim, in the text of the paper, 
that Erzberger's test is a necessary condition, whereas it is only 
sufficient, and quite restrictive at that. They do, however, use his 
results properly, which relegates their error to editing. Their approach, 
as optimal control approaches go, seems reasonable. At least they try to 
take up some of the slack in the problem by introducing side conditions, 
rather than ignoring it and plunging blindly ahead. 

As a whole, then, the literature is none too satisfying on this 
problem. The solutions are sparse and essentially ad hoc . None is of 
the type where the problem may be simply packaged and fed to a computer 
which could return an answer. All have to interact with their programs 
in a series of cut and try iterations until, through chance or artistry, 
the solution emerges. It is believed that the present work eliminates 
these objections. 

C. Allied Problems 

The most closely allied problem found in the literature is that 
of decoupling [7], [15], [22], in such a problem the designer is given a 
system and required to find a compensator which would decouple its transfer 
function. This is generally interpreted as meaning that each input should 
effect one output (or group of outputs) and no other input would effect 
that output (or group). If the transfer function matrix were square, 
we would wish to diagonalize it. 

The chief difficulty in turning this into a model following 
problem is the lack of a specific model. It is possible to cut and try 
models until a suitable one can be found, but this becomes tedious. 
Nonetheless, the two problems are closely enough allied that workers in 
either area can profit from papers concerning the other. Perhaps the 
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most interesting contrast between the two is that the 'decoupling' 
literature is almost exclusively theoretical in nature, while the 'model 
following' papers are devoted to computation and algorithms. This difference 
is perhaps most valuable to the workers in decoupling since many model 
following algorithms could be used to decouple whereas little of their theory 
can be used in a model following algorithm. 

The present work has had considerable success in solving decoupling 
problems. Since 'exact' solutions are given, the results will indeed 
decouple the given plant. The only difficulty lies in finding a suitable 
model, in Chapter V we will show how the work of Wonham and Morse [15] 
and [22], is connected to the present discussion. 

D. Contributions of this work 

The essential contributions of this work are as follows: 

1. The introduction of and development of equicontrollability 
(Chapter III). 

2. Development of an algorithmic solution to the model- 
following problem which is suitable for computer 
programming (Chapter IV). 

3. Resolution of some questions on the decoupling of 
multivariable systems (along with the application of 
the above algorithm to accomplish such decoupling) 

(Appendix C) . 

Attendant on these specific items are the techniques devised for the 
analysis of multivariable systems - particularly the use of non-minimal 
(equicontrollable) realizations, as noted in Chapters III and IV. 
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II. PROBLEM DEFINITION AND DISCUSSION 


In view of the preceding discussion, it becomes apparent that 
some better approach must exist for solving the model following problem. 

It should not require a performance index as such, and should yield simple 
answers to simple problems. It should be reducible to an algorithm 
suitable for digital computation with little or no human interaction from 
input to solution. Hopefully, it should not be iterative. The solution 
to be proposed satisfies all of these requirements. 

As has been pointed out, the definition of the problem essentially 
structures the solution. We suppose that we are given two linear constant 
dynamical systems in ithe differential equation form: 


PLANT: 


x = F x + G u 

nxn nxm 


y = H x 
mxn 


( 2 . 1 ) 


MODEL: 


6 = ® i + r u 

pxp pxm ’ 


H = A I 
mxp 


( 2 . 2 ) 


All capital letters are matrices. Note that the inputs and outputs 
of the plant and model have the same dimension. This is not an unreason- 
able assumption since it is hard to conceive of a problem having a plant 
with, say, three outputs attempting to track a model with four! This 
would be contradictory. The requirement that the number of inputs 
equal the number of outputs is more obtuse, but not without justification. 
In most aircraft problems, the number of control inputs is equal to the 
number of consciously controlled outputs. That is, if the pilot pushes 
the stick forward he has in mind a certain response in pitch. That 
this will, some integrations later, produce a change in altitude is not 
of immediate concern. Better, if the aircraft executes the proper 
motion in pitch, the other modes attendant on pitch will respond favorably. 
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Since the problem itself is linear, we might require that a 
linear controller be used. Moreover, since virtually all other methods 
involve state feedback, it is logical for us to do so also. This may be 
reduced to specifying that the controller be a linear function of both 
the states and external inputs. In the block diagram form shown in 
Fig. 2.1 we insist on the linear transfer function C(s) as the compensator. 



Fig. 2.1 

Form of Compensation 

The choice of a specific controller C(s) Must be based on 
some measure of goodness. We choose to make this measure in terms of 
acceptable modifications of the model instead of the minimization of some 
function. Although the advantages of this stratagem are legion, the 
clarity with which any performance degradation is displayed before the 
computation is performed is, in itself, sufficiently compelling to 
justify its use. More explicitly, we shall define the new model in Fig. 
2 . 2 : 


co 



Fig. 2.2 
Form of Model 

where a is prespecified by the designer, and r is to be as small 
as possible . 
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a is a clear measure of closeness. As a -* oo , the above transfer 
function approaches V(s) . But the penalty for making a large is 
to have large feedback gains. Hence, one may conceivably balance gain 
magnitude against accuracy in matching, which is an anticipated tradeoff. 

This discussion may be maneuvered into a definition of the 
problem, by way of listing the assumptions we shall now make. 

1. The plant and model are finite dimensional linear 
constant systems. 

2. Both the plant and model are completely controllable 
and observable. 

3. The plant and model have the same number of inputs and 
outputs. (All four are the same). 

4. The compensator should be a finite linear constant 
system, C(s) , taking u , the command inputs, and x , 
the plant state, into u . This assumes that the state is 
available. 

5. s = 0 is not a root of either mam det H(sl-F) 1 G or 
num det ACsI-G) - "^ . (Num det W(s) is defined below). 

6. Num det H(sl-F) has no roots in the right half plane. 

7. The designer will accept a revised model of the fonn: 

-l a r 

A(si-s>) r — ~ 

(s+a) 

where a is at his disposal, and r is to be as small as 
possible. 

The assumptions, except for 5 and 6, are direct consequences of the 
previous discussion. The fifth is best taken on faith until Chapter 4 
in which we shall use it as a guarantee that the algorithm will work. 

The symbol num det W(s) (where W(s) is a square transfer function 
matrix) represents the numerator of the determinant of W(s) . The 
numerator is extracted when the denominator equals the characteristic 
equation of W(s) . A more graphic description is: if W(s) = 

H(sl-F) "^G , then: 

num det W(s) = det W(s)-det (sI-F) (2.3) 
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For transfer functions of (2.1) or (2.2), num det is a polynomial in s 
of degree at most n-1 . At the end of Chapter 3 we shall prove the fact 
that, for all conformable K , (i.e. for any K of the proper 
dimensions) num det H(sI-F)“ 1 G = num det (sI-F-GK) -1 G . In other words, 
num det W(s) is invariant under state feedback. 

Since, in the scalar case, feedback is unable to move zeros, we 
will have to use a series compensator and pole-zero cancellation to 
achieve a perfect match if plant and model zeros differ. Clearly to 
cancel an RHP zero requires an unstable root, which would result in 
an unstable design. Hence the inclusion of Assumption 6, which is the 
analog of the above observation in the multivariable case. This 
assumption can be weakened, as we will see later. 

Defining the problem via a set of assumptions is perhaps the 
clearest approach in this case. The original statement is so vague that 
no straightforward definition is really feasible. We do gain something 
for our interpretation of the problem. The ability to give the designer 
a remarkably lucid picture of what he can expect by way of a solution 
is especially useful. We do not have to present a vague idea of how well 
the compensated plant will do — we can nearly describe its transfer 
function! Only the value of r is unknown. The remainder is specified 
by the given equations or is up to the designer's whim. To ask for 
more would border on the unreasonable. 
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III. EQUI CONTROLLABI LI TY 


A. Definition and Motivation 

... . .. . 

For a number of problems, the usual property of controllability 
does not provide sufficient information as to structure. A controllable 
system with m inputs may be controllable from only one input or require 
each and every one. Moreover, one cannot specifically give a canonical 
form for multi-input systems that is the same for all such systems. The 
best you can do is specify a procedure which will always yield some form, 
although many forms are possible. We shall consider here how we may 
trade one form of uncertainty for another to best serve our purposes. 

DEFINITION : The pair (F,G) has controllability index p 

if p is the smallest integer such that: 

C p = ^G,FG, . . . ,F P_1 gJ (3.1) 

has full rank n = dim F . 

DEFINITION : The pair (F,G) is equicontrollable if n = pm 

where p is the controllability index of (F,G) , F is nxn and G 
is nxm . 

Note that the definition of equicontrollability implies that 1 
m| n and that the first n = mp columns of the controllability matrix 

1 i 

m| n is read "m divides n" and means that there exists a k such that 
n = mk where k belongs to the same class as m and n , (e.g. if 
m and n are integers, k is an integer). 
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are independent. These conditions are enough to enable us to write such 
a system in a particularly useful canonical form — equicontrollable form. 
Let us consider a slightly more general case. 

PROPOSITION 1 : Given the pair of real matrices (F,G) , F nxn , 

G nxm and full rank, then there exists a similarity transformation T 
such that: 


TFT*" 1 = 


m 

'e 


n-m 

I - 


L x J 


and 



0 
e 
e 
x 

L x 


TG = 


e 

e 

9 

x 

x 


m 

'0 


L S J m 


n-m 


r m-r 

0 0 


0 

0 

I 

0 ' 


0 

0 

0 

I 


0 0 

X X 


0 

0 

0 

0 

0 

I 

Y 


X J r 


m-r 


m-r 

r 

m-r 


’ m-r 

r 


~0 

0 

r 

► 

e . 

0 

m-r 

0 • 

0 

r 

i 

0 

m-r 

_A 

I _ 

r 


( 3 . 2 ) 


( 3 . 3 ) 


where the X’s, Y, and A are unspecified entries determined by (F,G) ; 
and the parameter r is defined such that 


n = (p - l)m + r 


0 r < m 


( 3 . 4 ) 


1 All unmarked braces represent blocks of m rows or columns. 
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if and only if the last n columns of the controllability matrix; 


C = F G , F G , . . . , FG , G 


( 3 . 5 ) 


are linearly independent. 


NECESSITY : Suppose the matrix T exists as specified and 

define the following matrices: 


F* = TFT 


G* = TG 


-1 


let 


C = last n columns of C 

c, = [! , r i o,. p r 2 °* f . g »- g .] 

= last n columns of C,. 

* * 


then by inspection, we see that 


and hence that 


C = T - 1 C, 


C = T - 1 C, 


( 3 . 6 ) 

( 3 . 7 ) 

( 3 . 8 ) 

( 3 . 9 ) 

( 3 . 10 ) 

( 3 . 11 ) 


Moreover, C is nonsingular by construction. More specifically: 


C * = 


r m-r r m-r r 

i e e e e 


x i e 

x A i 

xxx 


e e 
e e 


i 0 

x x x • • • A i 


m-r 


m-r 


( 3 . 12 ) 


thus C is nonsingular as required. 
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SUFFICIENCY : The proof is constructive. Assume that C is 

nonsingular. First we define partitions necessary to describe C in 
terms of F and G . 


Let 


m-r r 


° = [ G 1 i g 2 ] 


(3.13) 


then 


and 


C = 


G = [ FG 2 i G l] 


m 


; A 2 = 

i Fk ' 2 °l 

A 


m 


_k-2_ 1 

F g 2 : 


F k " 3 G i 

^ 1 "V" * 

m 


m 


FG 


(3.14) 


m 


. , FG t FG_ ; G 

i -Li & , l 


G 

m 


— 3- 

r 

(3.15) 


Next we define the mxn matrix E from which the transformation is to 
be constructed. 

mm m r 

EC = [ I 9 . . . 9 9 ] m (3.16) 


E exists since C is nonsingular, and is simply the first m rows of 


-r-1 


We partition E so that: 


n 


E = 


L- E 2- 


(3.17) 


m-r 


We claim that the transf ormation T is: 


T = 


E 

EF 


' k— 2 
EF 


E F 


k-1 


m 

m 


m 


(3.18) 
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In the first place, T is nonsingular as we can see by using the last 
partition given for C : 


By the definition of 
to give: 


TC = 


k-2— 




EF G . . . 

EG 

EG 2 

k-1— 




EF G . . . 

t 

EFG 

EFG, 

• 

• 

2k-4- 

„ k-2- 

k- 

KK Q « • • 

EF G 

EF 

„ 2 k- 3— 

k-1— 

] 

E^F G • • • 

E F G 

E X F 

elements in 

this matrix 

can 


(3.19) 


k-1. 


2 


TC 


EF 


m 

I 

k-1- 


m 

0 • * • 


m 

9 


0 

0 


m 


m 


(3.20) 


2k- 4- 
EF G 

„ 2 k- 3— 

E F G 


E 1 F k_1 G 


0 

I 


m 


Hence T is nonsingular since C is nonsingular by hypothesis, 
we find the product: 

m-r r 


E 1 G 


0 

© 

r 

► 

E 2 G 


© 

© 

m-r 

• 





• 





E 1 F k “ 2 G 

- 

0 ' 

© 

r 

E F k_2 C 


I 

0 

m-r ' 






E F k-1 G 


A 

I 

r 

1 _ 


L- 

— 



Also 


(3.21) 
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Then as a special case we may write 


COROLLARY: Given the conditions in the previous proposition, 

and given that r = 0 (i.e. m| n) then there exists a transformation 
T such that: 



(3.23) 


iff (F , G) is equicontrollable. The above is equi controllable canonical 
f orm . 

The name "equicontrollable" is now easy to justify. Fig. 3.1 
shows a block diagram representative of such a system. Each input drives 
a string of m integrators, hence in some sense, the inputs control 
"equally”. The figure also shows (as does the canonical form) the 
similarity of this form and control canonical form. 

PROPOSITION 2: The form given in Proposition 1 is unique. 

- S /Si /\ 

That is, given two realizations of the same system: (F,G) and (F,G) , 
both in the canonical form specified in the statement of Proposition 1 , 

A A 

then F = F and G = G . 
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The proof consists of two lemmas. 

Lemma i : The recipe given in Proposition 1 yields the same 

result regardless of the coordinates of the original system. That is, 
given (SFS , SG) , the procedure of Proposition 1 yields the same 
realization, for all S . 


Proof of Lemma 1 : , In the sufficiency argument of Proposition 
1, we constructed the transformation T which produced the desired 
realization. If we define 

(F S , G S ) = (SFS -1 , SG) (3.24) 

and: 
also: 

and so forth, we find that: 
so that: 

s s 
T F 

and 


Lemma 2: If (F,G) is in the canonical form given in Proposition 

1 then the transformation T given in the proof of sufficiency equals I . 


Proof of Lemma 2; Consider C . After some tedious computation: 


5 s 4 [(F S ) k - Vi (F S ) k - Vi 


G S] 


(3.25) 


E S (? = [I 0 ... 0] 


(3.26) 


C S = SC , E S = ES* 1 , and T S = TS' 1 (3,27) 


, s -1 - Is— 1 ___-l 

.T ) = TS F ST = TFT 


(3.28a) 


T S G S = TS _1 G S = TG 


(3.28b) 


Q.E.D. 
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where the 


r 

m-r 

r m-r 

r 


r 

m-r 

r 


i 

! e 

i 

0 

e 

0 


0 

0 

0 

" 1 

e 

i 

• i 
j _ . 

1 

0 ; 
’ " “i 

0 

0 


0 

0 

0 

m-r J 

X 

! A 

i 

i ; 

0 

0 



• 





r 





* 



X 

X 

x i 

1 

i 

0 ! 

i 



» 



X 

X 

1 

x ; 

A 

i 

i ! 
— r 

t 

1 

I 





X 

X 

X 



• 

I 

0 

0 

r 

X 

X 

X 


• • 

• 

0 

I 

0 

m-r 

X 

X 

X 




X 

A 

I 

r 

's 

are 

unspecified entries. 

This implies that 




m 

m 

m 

m 

r 





E = 

[I 

0 

0 ... 

0 

0]m 





EF = 
« 

[9 

I 

0 ... 

0 

0]m 




EF 

• 

« 

k-2 

[0 

0 

0 ... 

I 

0]m 




E F 

k-i _ 

[0 

0 

0 ... 

0 

I]r 




(3.29) 


(3.30) 


hence T = I . 


Q.E.D. 


Proof of Proposition 2; Lemma 2 shows that the recipe in 
Proposition 1 yields T = I if the realization (F,G) is already in the 
given canonical form. But Lemma 1 shows that the procedure of Proposition 
1 always gives the same answer, regardless of the coordinates of the 
original system. Hence any two realizations in that canonical form 
must be the same. 

Q.E.D. 
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COROLLARY 1: Given two minimal equicontrollable systems of 

the same state dimension (H, F, G) and (A, $>,D , there exists a 
feedback matrix K such that: 

H(sl-F-GK) _1 G = A(sI-<J>) -3 T (3.31) 

A A 

if and only if H=A when the two systems are reduced to the respective 
equicontrollable canonical forms: (H, F, G) and (A, $, D . 


Sufficiency is easy. If the last m rows of F are [F ,...,F ] and 
of $ are J then the feedback matrix K is: 

K = [4, - F, - F IT" 1 (3.32) 

11 P P 

A A — 1 1 /\ A 

where T$T = 0 ; TF = r , AT = AT = A (= H by hypothesis). 


So 


T(F + GK)T 


A 

= <D 



(3.33) 


Hence, via the proposition, we are done. 

Necessity is similar. Suppose K exists as required. 
Then there exists a T such that: 


wr" 1 = F + GK 


Tr = g 


(3.34) 


~ -1 

AT = H 

But (H, F + GK, G) and (A,0, f) satisfy the hypothesis of the 

BHI 

proposition, hence T = I and A = H . sjjsj 

COROLLARY 2: In proposition 1, Y = 9 , i.e.: 


25 - 



TFT 


-1 


e e i © 

e e 0 i 

• 

e e 
e e 


x 

x 


X 

X 


e e 
e e 

I 9 
0 I 
X 0- 
X X 


N.B. 


(3.35) 


and 


Hence: 


but : 


so: 


Proof : 


Y = E 2 F k-1 T -1 


0 

0 

I 


tg 2 = 


’Q~ ' |m 
« 

0 |m 

Li. 


Y = E 2 F k “V 1 (TG 2 ) = E 2 F k_1 G 2 


k-1 

EF G 2 = 


0 


m-r 


Y = E 2 F k_1 G 2 = [9] m-r 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


(3.40) 

Q.E.D. 


As an interesting observation, if m| n and the system 
has been placed in the canonical form given in Proposition 1 but Y ^ 0 , 
then the transformation 
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r m-r r 





(3.41) 


will leave the form the same, only affecting the last m rows of the 
system matrix. Hence, the form is surely not unique. That is, other 
transformations than (3.18) will yield a similar form. (Note that if 
r = 0 or Y = 0 , then Q = 9). 

The utility of equicontrollable form is yet to be demonstrated. 
The reader could well interject that for n = mp to be true requires that 
a substantial structural constraint be placed on the system. To widen 
the class of minimal realizations which can be placed in equicontrollable 
form has been shown above to be impossible. On the other hand we could 
allow the addition of states to the system, which do not modify the 
transfer function, but do alter the structure so as to achieve equi- 
controllability . Such states will have to be controllable but not observ- 
able. (One gets into some semantic difficulty here. Since we shall be 
physically constructing such states — say in analog computer fashion — 
they are "measurable", but from the outputs specified , they are not 
observable.) The following theorem will demonstrate that we can always 
add such states. The method to be used will then be displayed. 

THEOREM 1 ; Given a system (F, G) with controllability 
index p, F nxn , and G nxm : 

If n < mp , there exists matrices A and B and a number 
N such that the N-dimensional system: 
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(F, G) 


(3.42) 


is equicontrollable iff N = m(p + k) for some integer k ^ 0 . 


SUFFICIENCY ; The proof will be essentially constructive. 
First the system is reduced to Luenberger form [13]. Fig. 3.2 shows 
a block diagram of the system, where we have renumbered states so as 
to draw the integrator strings in order of increasing length, left to 
right. Now we simply add integrators to the blocks so as to make them 
all the same length. 

The smallest number we may add to do the job would be 

J5, th 

\ (q -q.) where q. is the number of integrators in the i — block, 

Zu m 1 1 

i=l 

and q is the length of the largest such block. But; 


I 


(q -q,) = m 

m i 


m 


2 ’* 


= mq - n 
m 


(3.43) 


Moreover, if we added another k integrators to each block we would have 

m(q +k)-n states added, for a total of m(q +k) states after all have 
m m 

been added. Hence, if we can show that p = q^ , we are done. By the 
definition of p , 


C = [G, FG 
P 


F P-1 G] 


(3.44) 


will have n independent columns. Recall that in constructing the 

th 

Luenberger canonical form, the length of the i — block was determined 

k- 1 

by how many vectors of the form F q^ were independent of vectors 

previously chosen from the controllability matrix. However, by the 

definition of C , k £ p regardless of how the vectors were chosen, 

i.e. q < p , but suppose q < p . Then k < p which says that no 
m ' m 

block has length p so that 
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q -i 

C = [G, FG, F “ G] 

q m 


(3.45) 


has full rank! But this contradicts the definition of p , hence 

q = p . 
m 

NECESSITY ; We suppose that A and B exist, and that the control- 
lability index of (F,G) is p . If we show that p ^ p , then 
N mp . By appealing to Proposition 1, m| N so N = m(p+k) for k 
an integer 0 . 

Suppose p < p . Consider the partial controllability 

matrix: 


= 

P 


G FG F G 

e AG AFG+BAG 


A 

.. F P G 


X 


(3.46) 


C/>, is NxN and full rank. Therefore, 

P 


[G, FG, ...,F P G] 


is full rank. But this contradicts the definition of p , hence 
0 p , and N = m(p+k) as previously argued. 


To demonstrate the nature of the sufficiency argument (and 
lead into an algorithm for implementing it) consider this example. 


F = 


r_ 




— 






f 1 — - 





___ cr * m 





f ' _ 




g l 

fr • « 




2 

- f ' 



— R 1 






3 



B 3 



i 

0 0 

0 

0 

0 

, G = 

0 

0 

0 

0 

1 0 

0 

0 

0 


0 

0 

0 

0 

0 0 

1 

0 

0 


0 

0 

0 










—* 


(3.47) 


The system may be represented by the block diagram in Fig. 3.3 
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f.|x+gju f^x+g^u 




f *x + g' u 





ORIGINAL 

SYSTEM 



/ 

\ 


I ADDED 
/ DYNAMICS 


J 


Fig. 3.3 
Block diagram of 
equicontrollable form 


The resulting time domain matrices (after dynamics are added) are: 







— 

* 0 
1 

0 

o” 







_f 

t 

2 - 


— 

1 

i o 
1 

0 

0 









_f ' 


1 

1 0 

0 

0 









3 


1 





G 


1 

0 

0 

0 

0 

0 

! ° 

0 

0 





0 

1 

0 

0 

0 

0 

1 

1 0 

0 

0 











1 



, G = 




0 

0 

0 

1 

0 

0 

i 0 

0 

0 





— 

— 

— 


— 

- 

1 — 

- ' 

— 


— 

— 

— 

0 

0 

1 

0 

0 

0 

» 0 
1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

1 

0 

I 0 
1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 1 

0 

0 


0 

0 

0 







1 




— 


— 


(3.48) 


Now the system may be put into equicontrollable form by renumbering 

the states. (Recall p = q hence, since n = q , n = p ). 

m mm 
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Although not explicitly stated, the method of adding states 
used in Theorem 1 works equally as well if, instead of each added state 
being appended as — , it is appended as — . A glance through the 

S S 401 

proof of sufficiency will confirm this. Since the resulting system is 
equicontrollable, we know via Proposition 2 how to reduce it to equi- 
controllable form. 

In summary, we have shown: 

1. That n = pm is the key to equicontrollable form. 

2. That there is some measure of uniqueness to such 
a form. 

3. That we can find a non-minimal but nonetheless 
equicontrollable realization for any transfer 
function. 


B. 


Properties 


PROPERTY 1 : If the system (H, F, G) has been converted 

to the equicontrallable form: 


H= 


= [H 1 ,H 2> . . . ,H ] , F = 


“e i e ... e" 


e 

e e i ... e 


e 

• 


• 

• 

, G = 

• 

• 


• 

e e e ... i 


e 

F„ F n F„. . . F 


i 

12 3 p 




(3.49) 


Then W(s) = H(sI-F) -1 G = 


"f v 1 ' 1 

s K I - \ F.s 
1 

Li=l j 

L i=l J 


-1 


The last column of (sI-F) 1 is easy to compute. To wit: 
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W(s) = 


* * * ,H p 


I “ 

si 

lx 


.1 XU 
s' i - y f ± s 

i=l 


i-1 


1-1 





p ^ i-1 

sr 

\ H.s 


s i - \ s F. 


Z 1 


Z 1 


_i=l j 


L i=l J 


-1 


(3.50) 


PROPERTY 2: Num det W(s) = det 


i=l 


. s 


i-1 


= num det 

r £ i-ii 


P i-1 

V H.s 


s H I - \ F.s 


Z 1 


/ l 


L i=l J 


i=l 


= num < 


f £ 

Idet y 

„ i-l 
H.s 

l 

det 

S P I - 

1 

„ i"l 

F. s 
l 

-l 

l 




i= 

i 

> 

det 

P 

I H ‘ 

i-1 

s 


\ 




i=l 






det 

s p I - 

Z i-1 

\ F.s 

A 1 

J 

► 



(3.51) 


l (3.52) 


i=J 


(3.53) 




= det \ H.s 
i=l 


i-1 


(3.54) 


The last step follows from the definition of "num det", 
equation (2.3) and the observation that: 


X„(s) = det (sI-F) = det 

r 


s P ! - ^ - 1 " 1 ' 


X- x 

1 8 F x 


i=l 


(3.55) 
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PROPERTY 3: For any conformable matrix K , 


num det H(sI-F) _1 G = num det H(sI-F-GK)“ 1 G . 


(3.56) 


PROOF: (cf: Morgan [14]). Without loss of generality, we 

may assume that: 


F = 


F , 

F, „ 


"e 

11 ' 
i 

v 

12 

, G = 


F_ • 

F 


i 

[_ 21 i 

22 J 




H =[ H l| \ 


(3.57) 


Let: num det H(s-F) ^G = num det W (s) and T = 

F 


-i : < s i-f 11 )- 1 p 12 ' 1 


det W ( s ) = det HTT* 1 (sI-F)” 1 G 
J? 


(3.58) 


= det 


[ H lj H l <8I - f 'll ) ' lF 12 +H 2] 


' 

sI-F . ! 

-F ~1 


11 J 

12 


L 


1 

-F_ ! 

Sl-F_„ 


L 21 1 

22 J 


: <si - p u rlF i2' 

- 1 

e 

1 

j 


i 

> 


\ (3.59) 


= det 


[ H i| H i 


(sI - F U> F 12 +H 2 


sI-F 


11 


-F 


12 


(SI_F 22 )_F 12 (SI F ll ) F 12 


-1 


9 


(3.60) 
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(3.61) 


det [ H l ( s 1 - F 1 l ) ~W H 2 ] 


det f 


Sl-F 22 -F 21 ( sI-F 21^ 


' X F 1 
12 ] 


Observe that: 


1. detf^Ul-F^F^H,,] = 


det 

si - f h : 
r 1 

_” F 12 


H i : 

H 2. 

det 

(sI_F ll ) 



Therefore: 


det 


det W(s) 
F 


11 ) . det 

(sI-F 

sI-F . ' 


11 1 

12 

H_ ' 


L 1 

2 J 


22~ F 21 (SI-F 11 ) F 12^ 


X F (s) 


so that: 


num det W ( s) = det 
F 


sI-F 


11 


H. 


12 


H 


2 -» 


But if we compute num det W (s) in the same fashion, 

F+GK 


F+GK 


sI-F 1 -F 
11 • 12 


= det 


det [ Hl !H 2 J 

-F - 
L s 21 



sI-F. . ! 

-F, „ 

11 1 
1 
1 

12 

1 

H 1 

H„ 

L 1 

2 J 


-1 


22 2 J 


num det W_(s) 
F 


Q.E.I 


(3.62) 

(3.63) 

(3.64) 

(3.65) 

" 9 " 

I 

(3.66) 

(3.67) 

(3.68) 


35 - 



Note that property 3 is not only a property of equi controllable 
systems, but of any system. We include it here since it fits in well 
with the properties of equicontrollable systems, making a well rounded 
picture. 

PROPERTY 4; Suppose (H,F,G) is equicontrollable and 
completely observable, and (C,A,B) is equicontrollable. Further, F 
is nxn, A is NxN, N ^ n , and both have m inputs and r outputs. 
If 


C(sI-A) _1 B = H(sI-F) _1 G 


(3.69) 


then there exist T,X,Y such that: 

[H ;9]T = C 
G 


-1 


.0 


f ; e 

l~ 


x ; y 


= B 


T .= A 


(3.70) 


Proof: Observe that since the transfer functions are equal, 


k k 

HF G = CA B for any k 


(3.71) 


Let p and q be the controllability indices of (F,G,H) and (A,B,C) 
respectively. Then we claim that: 


T = 


G FG 


F P-1 G ! F P G 


f’^g 


^ ' 
I 
I 
I 


m(p-q) 


[ 


B AB 




H ' 1 


(3.72) 


It is easy to check the first two equations above. To wit: 
£ H ; ©Jt = ^HG HFG .... HF q_1 Gj |b AB . . . A Q_ 1 bJ 


(3.73) 


but from our initial observation, we may write: 
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Finally we must check the last equation. 



ASIDE: From the above: 
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( 3 . 79 ) 


A q B 


-[■ 


B AB 


i'-S] 


L U 


Then; 


c 


C 


‘V 

CA 


CA 

r . „ . q- 1 1 

A. 


AB = 

■- 

B AB ... A^ B 

.2 

• 

• 

A 





L q J 


which implies that: 


- - 


■» — 



H 


H 


’V 

HF 


HF 

r q-1 "1 

A_ 

• 

• 

• 

, n-1 
HF 

F G = 

• 

• 

• 

t n-1 
HF 

G FG ... F* G J 

.2 

« 

A 

*— - 




L q J 


( 3 . 80 ) 


( 3 . 81 ) 


But Q is full rank (since (H,F,G) is completely observable). Hence 
we conclude that; 


q 

F h G 


-[■ 


G FG 



A 2 


A 

q 


( 3 . 82 ) 


Now to return to the original train of thought. We may consider the first 
n rows of the matrix TAT 1 as follows. 
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j^G FG ... F Q_ 1 gJ 


Hence: 


0 0 • • • 0 
10 0 A 2 

0 1 0 A„ 


0 0 • • • I A 


_ -1 


= j^FG F 2 G ... F q_1 G F q Gj 


- q-1 

G FG ... F 4 G 


0 ; i 


q- 1 — | 1 

G FG ... F 4 G 


(3.83) 


= fJg F 2 G ... F q-1 Gj 

= f '[ I ni 9 ] 

= [f; 0 ] 


G FG ... F q_1 G 


TAT 1 = 


rF 1 0l 

I 

- - 
I 

x 1 y 


for some X and Y . 


(3.84) 

(3.85) 

(3.86) 

(3.87) 
Q . E . D . 
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IV. THE ALGORITHM 


A. Introduction to the Algorithm — The Scalar Case 


The basic difficulty in presenting the algorithm is in the 
large number of steps required. Any explanation tends to lose continuity 
when attempting to clarify each individual step of such a long 
procedure. To attempt to alleviate this problem, we shall try to 
proceed in a number of progressively more difficult stages, each 
cutting deeper into the wealth of detail involved. Extensive use will 
be made of simultaneous presentation of steps in sundry forms (block 
diagram, time domain equations, transfer functions, and text description) 
so that the mechanisms involved will be as clear as possible. 

To begin the discussion let us consider a scalar problem in 
transfer function terms. The model and the plant might, for example, 
be given by: 


PLANT 


0 


(s+1) (s-2) 
(s+3) (s+4) (s+5) 


MODEL 


0 


s+3 

(s+6) (s+7) 


(4.1) 


A moment's reflection will reveal that no realizable 
compensator can be added to the plant which will cause it to behave the 
same as the model. Further, if we are to be able to find a stable 
compensator, we shall have to require a change in the model, namely it 
must also have a zero at s = 2 . These observations may be summed up 
by this modified problem. 
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PLANT 1 : 


(s+1) (s-2) 
(s+3) (s+4) (s+5) 


MODEL 


s +3 

(s+6) (s+7) 


g(s-2) 

(s+a) 


(4.2) 


a is arbitrary, but is prespecified by the designer. If this model's 
response were not acceptable, there would have to be negotiation with 
the designer until a suitable model were found. But we suppose here 
that this model was approved. Then the way to a solution is clear. 
First, we find the feedback which matches the dynamics (poles) of the 
two systems, yielding: 


PLANT 


2 


(s+1) (s-2) 
(s+6) (s+7) (s+a) 


MODEL 


2 _ 


s+3 a(s-2) 

(s+6) (s+7) * (s+a) 


(4.3) 


Finally, we apply an input series compensator to the plant, which 
completes the solution. 


PLANT 3 : 


(s+1) (s-2) 
(s+6) (s+7) (s+a) 


as(s+3) 

(s+1) 


MODEL 


3 


(s+3) a(s-2) 

(s+6) (s+7) ’ (s+a) 


(4.4) 


/ 

The transfer functions are now the same, so we are done. Some comments 
are in order at this point. 

1. Although this was a specific example, the procedure 
is quite general. Even if the plant had fewer states 
than the model, the same method would be used. 

2. There is no doubt about the minimality of the solution, 
since the manipulations are quite transparent. 

3. The need for adding RHP zeros of the plant to the model 
is also obvious, lest the compensator be unstable. 

Let us do the same problem in the time domain. The format will be to 
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place the associated block diagrams, time domain equations, etc., in 
an accompanying figure. Thus the steps may be more easily visualized 
during the discussion. 

We begin with the original problem (Fig. 4.1). Since the 
systems were assumed to be controllable, we have chosen to write them 
in control canonical form. (Ordinarily the problem would be given in 
differential equation form, so it would first have to be reduced to 
such a form.) Exactly as before, the first step is a modification of 
the model. 




Fig. 4.1 

Scalar Example: 'O' Coordinates 
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Fig. 4.2 shows a method of appending the additional state to 
the model. This particular way of introducing them is used since only 
input changes are deemed reasonable (or realizable) . Although such a 
restriction is not strictly necessary in the model, which is a 
mathematical fiction, it is essential in dealing with the plant. The 
outputs represent (usually) physical actions, such as pitch or roll 
of an aircraft, which must be the outputs. This observation is basic 
to the problem but has occasionally been overlooked, such as in [2]. 
Generally, this step would also include another transformation to 
control canonical form, but it was not shown here for the sake of 
better displaying the details of the state addition. 


















At this stage (Fig. 4.3), three changes have been made. 

First, the model state has been redefined so as to place it in 
control canonical form. Second, the plant has had feedback (per 
fig. 4.4) applied, viz: 

u = u + Kx (4.5) 

Third, a series compensator has been added. In the figure, the 
feedback is evidenced by the third row of the F matrix, to the 
left of the dashed line. The series compensator has state z and 
transfer function a(s+3)/(s+l) . The result has been to give the 
plant the same transfer function as the model (after revision). 

The compensator is stable, as is the compensated system as a 
whole. Moreover, the compensator is minimal, whereas the resultant 
system is not, due to the cancellation of the (s+1) factor. 

Fig. 4.4 shows the series and parallel compensation as applied 
to the original plant. Note that, in this case, there is no 
feedback from the states of the plant to the input of the series 
compensator. This is not a general property of the method. 

The changes in the model are neither unexpected or 
unreasonable. Were it not for the RHP zero of the original plant, 
the changes required would have been especially acceptable. 


44 



2 

PLANT : 




o 

o 

o 

1 


~o ~ 


• 

A 

X 


o 

o 

o 

A~ 

X 


0 


• 


1 

— 

+ 


C> 

II 

i 

i 

i 

A 

Z 


-42 -(42+13a) -(a+13) ! 1 

A 

Z 


a 




0 0 0 !-l 


2a 





A 


\ 

y 

= 

-2 -1 1 '0 



A A x 

= H 




z 


z 




- 




W(s) = 


(s+1) (s-2) 
(s+6) (s+7) (s+a) 


g(s+3) 

(s+1) 



MODEL 


2 


unchanged 


V(s) 


(s+3) 

(s+6) (s+7 ) 


a(s-2) 

(s+a) 


Fig. 

Scalar example: 


4.3 

’2' Coordinates 
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K = [ -42a+60 j - (42+1300+47! - (ot+13)+12 ] 

A-i 

u 

z = [-l]z + [2*9] — 

* x_ 

u “ 

u = [ l]z + [ct;K] •'* 

x _ 

Fig. 4.4 

Scalar Example; Compensator Structure 

* 

In summary then, the steps involved in a scalar problem may 
be given as; 

1. Get both the plant and the model into control 
canonical (equicontrollable) form. 

2. By adding states to the plant and/or the model, correct 
unequal numbers of states in the plant and model. 

3. Adjust pole locations by feedback. 

4. Cancel and/or replace zeros in the plant as required, 
by series compensation. 

5. Transform back to original coordinates and write the 
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compensator as: 


z = F z + G 
c c 



u = H z + J 
c c 



(4.6) 


thus completing the design. 

Although the detail differences between the scalar and multi- 
variable problems are legion, the basic steps are very similar, and the 
result is expressed in the same form. 

t 

B. The Multivariable Case — Overview 

The problem to be faced when looking at the multivariable case 
is the lack of a good canonical form for all possible systems. Since 
no one is suitable, the approach we shall use is to increase state 
dimension to simplify structure. More specifically we shall add states 
to any system which is not equicontrollable so as to make it equi- 
controllable. In chapter 3, theorem 1, we showed this to be always 
possible and subsequently described the procedure that we shall adopt 
to accomplish it. We lose simplicity in the sense of state dimension 
and we introduce some non- uniqueness, but we are more than repaid in 
structure by being able to write every matrix in a canonical form having 
all mxm partitions. Moreover, the sundry properties of the equicontrollable 
form may be invoked. 

Even with this approach in mind, the path is not at all clear. 

A number of directions present themselves which seem at first to quickly 
solve the problem, but which actually create more difficulties than they 
remove. In this section we shall consider one of these, which, although 
it does not succeed, provides useful insights into the true nature of 
the problem and how we shall finally go about solving it. r 

To demonstrate the approach let us consider a simple multi- 
variable problem (Fig. 4.5). Both systems have been chosen to be equi- 
controllable and of equal state dimension. Thus the detailed (but 
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straightforward) step of achieving this form has been temporarily 
ignored. We intend to show that simply operating by analogy with the 
scalar case is not sufficient, and to do so would lead us astray, even 
in this simplified case. 

For each system the transfer function's form is derived via 
Property 1, Chapter III. 

The first step is a natural analog to what was done in the 
scalar case, but with a slight variation. If we wished to try a 
direct application of the scalar method, we would find a K such that: 

W s) = V s) < 4 - 7) 

This has two serious drawbacks: 

1. The feedback is not unique (as it was in the scalar case), 
so choosing it now would remove a degree of freedom in 
the solution. Thus it might be better to wait until we 
have better grounds for the choice. 

2. If we did choose K now, the next step could be impossible. 
Mainly since we will later add series compensation to the 
input , and we will need to have 

r P • 1 ‘i - 1 

s P I - ^ (F i + K i )s 1 
i=l 

and the compensator's transfer function commute. This 
is not always possible in the multivariable case, 
consistent with the requirement of prescribed character- 
istic equation. It was in the scalar case — a fact which 
was not explicitly stated. 
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In view of observation number 2 above, we can choose a K which will 


allow the necessary commutation. In particular choose K = -F ! 

i i 

The model is left unchanged, and the plant is still equicontrollable 

(a property invariant under feedback). Fig. 4.6 displays the result. 

Note that had H. = A.V i , we could have let K = 0. - F and been 
ii ill 

done at the first step! The chances of this occurring are so slim that 
we ignore it. Now, at least in principle, the way is clear. We build a 
realization of 


C(s) 


1 

Ip 

H- 

1 

-1 

^ i-1 

\ H. S 


\ A.s 

A 1 


Z 1 

Vi 

1 


i=l 


(4.8) 


if possible (we fail if num det W(s) has a pole in the RHP, for 
example, since then the compensator is unstable, recalling that 


num det W(s) 


divides 


det 



H 1 s x " 1 ). 


Then put that compensator in 


series with the plant to yield 

P 


W(s)C(s) = 


S 

i=l 


H i S 


i-1 


V H.s 1 - 1 

>4 1 


-1 


V 


i-1 


(4.9) 


i=J 


1*1 


i-1 


s P I| 


-1 


(4.10) 


i=l 


Now it would seem that we need only apply the feedback K^= and we 
would be done. Unfortunately, this argument is specious. It is not 
true that after application of the series compensator C(s) , we may 


apply such a feedback. 


In particular, the cancellation of 



i=l 


i-1 


in eq. (4.9) above does not occur internally, hence the state has increased 
in size and K cannot be as claimed. 


50 - 



PLANT 1 : 



MODEL 1 : 


unchanged 



Fig. 4.6 

Multivariable Example; ’1' Coordinates 

To demonstrate that this is the case, suppose p = 2 . Then 
say that C(s) has the (possibly nonminimal) realization: 

9 I 

[Bi B 2 ] , 

_ A 1 A 2 

(C(s) may require more or less states than shown but the principle 
will remain). Then the series compensated plant will have the form: 
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Now we are somewhat at sea. Is the system above equi controllable? 

Is it for some difficult realization of C(s) ? Even if it is, how 
do we find a feedback to match it to the given model? These questions 
do not seem to have general satisfactory answers at this point. 

On the other hand, the general approach is salvageable. We 
need only guarantee that a suitable feedback will exist to complete 
the compensation. This, in fact, can be done. In the next section we 
shall give such a construction. Briefly, it will proceed as follows: 

1. Place both model and system in equicontrollable 
form with the same state dimension, even at the cost 
of adding states. 

2. Apply a loop of feedback to the plant to place all of 
its eigenvalues at the origin. (In particular making the 
last m rows of its "f" matrix zeros.) 

3. Synthesize a series compensator along the lines of the 
method given here, but modified in such a way that indeed 
feedback alone will be able to complete the job. 

4. Find the aforementioned feedback, K . 

5. Untangle the compensator and return it to the original 
coordinates of the plant. 

6. Find a minimal realization of the compensator — considered 
as a linear system taking (x,u) into u . 

The next section will detail the algorithm. 


C. The General Multivariable Case — Detailed Description 

Fig. 4.7 depicts the various operations which we propose to 
perform on the plant. This is essentially a pictorial version of the 
outline just given, but containing substantial detail. Each step is 
numbered along the left margin. These numbers will correspond with later 
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numbering of the plant in various coordinate systems. Again, each of 
these coordinate systems will be considered in detail in subsequent 
figures. The transformations indicated on the right show how the' states 
of one version are derived from the previous one. In each case the 
block(s) within the dashed lines are derived through the indicated 
transformation from the previous diagram. The model will follow a 
similar history as we proceed, but it will not receive the various 
rounds of feedback. Fig. 4.8 gives the details for the model. 

Fig. 4.9 shows the starting place for the plant and model. 

Note that the figure gives the block diagram, a state description, 
and the transfer function for each. This pattern will be continued in 
the remaining figures. In general, the plant and model are of 
different dimension. Neither is equi controllable, and their controllability 
indices are usually different. To alleviate this uncertainty, we shall 
first add states to each to make them equicontrollable. (Alternatively 
one could say that we will find the smallest equicontrollable realization 
of each.) Then we shall increase the number of states in the smaller so 
that they both have the same state dimension. This addition will be 
made so as to create a new equicontrollable state realization of the 
appropriate state dimension. The reduction to equicontrollable form is 
accomplished in discrete steps. 

First, we add states with roots at the prespecified location 

cc . Fig. 4.10 shows the result of applying this procedure to the 

plant, as PLANT 1 . As the block diagram shows, the states z are 

totally isolated from the output. Hence they do not appear in the 

transfer function. J and J simply define the method of 

J. -L 

connection of those states, as per the recipe given in the proof of 
Theorem 1, Chapter III. In the example after that theorem, we may note 
that 



"0 

0 

1 

0 

0 

o' 


'0 

0 

o' 

J 12 = 

0 

0 

0 

0 

1 

0 

J ll = 

0 

0 

0 


_0 

0 

0 

0 

0 

0 . 


.1 

0 

0. 


2 

and a = 0 . PLANT is simply a redefinition of the previous coordinates 
where x and z are combined into the single vector x . 
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1 


2 


3 


4 


5 



Fig. 4.8 

Algorithm: Operations on the Model 
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% - + rw 

V = At 

V (s) = A(SI - $) _1 r 


Fig. 4.9 

Algorithm: plant and Model in 

'O' Coordinates 
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F5T + Gu 


x 


Hx 



Fig. 4.10 

Algorithm : Plant in '1' and '2 

Coordinates 
















Now the reduction to equicontrollable form, via the similarity 
transformation T , is shown in Fig. 4.1. In our notation, a subscript 
"*" indicates equicontrollable form, and will be referred to as "star" 
coordinates. 

The feedback KT sets the last m rows of F to 9 . We 
choose to define the feedback as KT^ for notational convenience later, 
when we will unravel the transformations to give the compensator in the 
original coordinates of the plant. The result of the feedback is that: 


F * + G*KT i = 


e i e e 

eel e 

• • 

• • 

0 0 0 I 

0 0 0 ... 0 


(4.14) 


Fig. 4.11 shows the external results of the feedback. The transfer 
function is reminiscent of the one found in the previous section (4-B) . 
From this point on it would be well to keep in mind the lessons learned 
in that development. 

Figs. 4.12 and 4.13 give the operations on the model so as to 
bring it to a set of coordinates comparable to the "bar- star" represent- 
ation of the plant in Fig. 4.11, but without the feedback. No special 
indication is made as to whether the plant or model had a larger state 
dimension at the outset. Whichever it was, it no longer matters. 

The discussion to follow is somewhat complicated, both in 
notation and in logic. Recall that if = A* , we could complete the 
problem at this point by feedback alone. We desire to find a way to 
add states to the plant and model so that, when again reduced to equi- 
controllable form, the output matrices will indeed be the same. More- 
over, we wish to do this consistent with the various restrictions which 
we placed on the original problem statement. Let us first look at the 
type of compensation we will add to the plant (series, in fact) and 
then go back and show how this answer works. Fig. 4.14 shows the 
application of the series compensator. J and J ^ are defined, 
and simply show how we shall wire in the additional states. We shall 

define k matrices A. and k scalars A.. by the boxed equation. 

i i 
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3 

MODEL : 



Fig. 4.13 

Algorithm : Model in '3' Coordinates 


61 - 






















PLANT 5 : 



Fig. 4.15 

Algorithm: Plant in '5' Coordinates 
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(Note that the right hand side of the equation is reminiscent of the 

C(s) given in section 4-B) . We wish to use the smallest integers r 

and k for which that equation is satisfied for some A. and X. . 

11 

The value of k is increased and zero A. are added if the numerator 

1 

order is too low. 

Fig. 4.15 completes the operations on the plant by applying 

a final round of feedback, KT . T represents the similarity 

A £ 

transformation used to transform the plant in Fig. 4.14 to the equi- 
controllable form in Fig. 4. 15. T is given by 

a 


T 2 = 


A 1 A 2 • • * A k 9 


e e 

e 0 


0 


0 A. A ... A 0 ... 0 
1 ^ k 


A 1 a 2 


(4.15) 


the 

matrices 

in Fig. 

4. 

14 

specifically as follows. 



0 

i 

0 • • • 

0 

• 

1 

0 

0 

. . . 

0 


0 



0 

0 

I 

0 

1 

1 

0 

0 


0 


0 




• 


• 

I 

• 



* 








. 

I 

. 



. 


• 

/N 

F = 


_0 _ 

0_ 

_0 J *_ 

_0 

I 

-L 

I_ 

0 


0 

o > 

II 

A k 



0 

0 

9 • • • 

0 

1 

0 

I 

0 • 

• • 0 


A, n 





1 




; 


k-1 



0 

0 

0 

0 

l 

l 

0 

0 

0 

I 


A o 





• • • 


i 

i 



• 







0 

0 

0 

0 

0 

0 

0 


_ A 1 

/\ 


r 




i 






H = 


1 H -, 

H. 

• * * 

H 

1 

0 

0 

0‘ 

- 0 



L 1 

2 



p! 




_ 



(4.16) 
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From which we may check that T„ is indeed the correct transformation. 

Note that T will be non-singular if and only if A is 
£ 1 

non-singular. That our assumptions are sufficient to guarantee the 
invertability of will be shown in Chapter V. 

We now claim that the compensated plant has the correct 
transfer function. To prove that this is the case, consider the 
further expansion of the model shown in Fig. 4.16. The transfer function 
has not been affected, but it has now as many states as the plant in its 
'5' form. Further, the state transformation: 


*1 

0 


*2 


0 


^2 


\-l 1 


0 

0 


S 2 - 


0 

I 

0 


0 

0 


V: 


0 

0 


0 

0 


0 

0 


I 

’0" 


-tH(4.17) 

i : © 


0 


I . 0 


will place that model into equicontrollable form. Now consider the 
following digression, which will show that: 


where: 


p+k-1 

I 


i=l 


H.s 

l 


i-1 


r 

a 


(s+a) r 


p+k-i ^ i _ 1 

s 


I V 


i=l 


/\ /\ A /N “J 

H * = HT 2 = L H 1’ H 2” ' * ,H p+k-lJ 

A * = ^2 = L A l ,A 2”*'’ A p+k-lJ 


(4.18) 


(4.19) 
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By rearranging equation (*) in Fig. 4.14 we obtain: 

i V 1 ' 1 - y v~ 1 


P . 1 k 
(s-*a) r ^ H i sl ~ 

1^1 i=l 


i-l r . i-1 ^ , i-1 


V 


i=J 


i=l 


equivalently: 


Cs-ta) r [H 1 ,H 2 ,...,Bl 


y\ 


i-1 


i=l 


IV 

Tv 


s \ A. s 
i=l 


i-1 


p-i V a i_1 
s > A . s 


Tv 

i=l 


= a 




IV 

2 v 

i=l 

k 

s X. s" 

Z 1 

i=l 


i-1 


i-1 


P-1 V' , -i-1 
i=l 


* 2 v 


(s+oj) iHg > • • • >^p 


= a 


[w-’Ap] 


A 1 A 2 A 3 **• \ 9 *** 
e A x a 2 a 3 ... A k e. 


_e e a x a 2 


\ *2 ^3 * * * \ 9 

9 *1 *2 V** \ Q ' 


9 

.9 


I 

si 


’ k+p-2 
s H I 


0 

0 


9 9 


\ ^2 




I 

si 


s k+p " 2 i 


(4.20) 


(4.21) 


(4.22) 
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Recalling that H = [ H| 0 ] and A = [a|Q] and the definitions of 
and we may then write: 


(s+a) r HT„ 


“ — 



I 


i 

si 

r A - 

si 


= a as 2 


. s p+lt - l I_ 


_s p+k_1 I_ 


(4.23) 


or: 


(s4a) r H. 


I 


’i 

si 

r /\ 

si 

_s p+k - 1 I_ 

= a 

_i p+k-1 !_ 


(4.24) 


p+k-1 . p+k-1 . , 

(s+a) \ h^s = a ^ a^s 

i=l i=I 


(4.25) 


But now note that: 


— 1 A A 

S 2 ** S 2 = 


9 I 
0 © 


0 ... 0 
I 0 


*1 *2 *3 Vk-1 


(4.26) 


By computing the above in a similar vein to the previous development, 
we can show that: 


p+k- 1 


I 


i-l 


i=l 


s P I - V o.s 1-1 


i=j 


r k 


I* 


i-l 


Li=l 


- s P+k_1 I (4.27) 


But this implies that = K. I So we may conclude that: 

ii 
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(4.28) 


A A _ l /\ 

H^CsI-F . ) G . 




= HT^sI-®^) 


-1 


— K^-h^r* 


(s+a) 
r 


— — a(si-<d) 1 r : 


(s+a) 


(4.29) 


(4.30) 

(4.31) 


Hence the algorithm gives the desired result. 

Unfortunately, we have two more tedious, but straightforward, 
steps to perform. First we must unravel the compensator and reduce it 
to original coordinates, and second we must find a minimal realization 
of that compensator. Both of these tasks are essentially bookkeeping and 
use standard techniques. The first may be accomplished by reversing the 
sundry steps shown in Fig. 4.7. The second can be achieved by using 
the technique given by Kalman [9], operating on the compensator shown in 
Fig. 4.17. 

The essence of Kalman's method is as follows. Suppose the 
given system were described by: 

x = Ax + Bu y = Cx (4.32) 


First check the rank of the controllability matrix: 


C = [B,AB, . . . ,A n_1 B] 


(4.33) 


If rank C < n , find the transformation T such that: 


T 1 C = 


i. 

Y 

» 

L e : 

0 J 


(4.34) 


Then: 


-1 

‘ B 

-1 

’ A • X " 

I 

l _ 

T B = 

“ ” “ * 

, T AT = 



_ 0 . 


_ e • y _ 


, cr = [c I c x ] 


(4.35) 


Therefore: 


T\( Q T A \ 


B(sI-A) -1 C 


(4.36) 
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Controller In Detail 



















One then replaces the system (A,B,C) with the smaller dimension system 
(A,B,C) . The same operation may be repeated in the dual, observability 
case. The final result will be a minimal (i.e. completely controllable 
and observable) system. The approach has no formal difficulty, but finding 
T accurately is no mean numerical analysis task, for even moderate 
state dimension. 

D. Conclusion 

By a somewhat circuitous path we have succeeded in solving, 
algorithmically, the model-following problem as posed in Chapter II. 
Although the process is involved, the operations at each step are 
quite straightforward. The resultant compensated plant has the proper 
transfer function, although it usually represents a non-minimal realization 
of same. Our solution differs markedly from the quadratic loss, but 
bears a great deal in common with, the classical approach. In a sense 
it is a hybrid, using the best features of each. One may visualize 
the process as finding feedback laws which take full advantage of existing 
plant structure so as to minimize the dynamics required in series. The 
classical approach insists on "unity" feedback, thereby obviating any 
potential advantage. 

Even when programmed in an essentially brute-force manner, 
the algorithm is capable of generating the compensator in about the 
time required for one pass through a quadratic loss program using the 
Q-R algorithm (a very fast technique). Considering that the quadratic 
loss method has only begun after one pass, our algorithm is especially 
attractive. 

Perhaps more significant is the potential of using non- 
minimal realizations (particularly equicontrollable or "equi-observable" 
ones) for the analysis of multivariable systems. Our previous claims 
as to the utility of such techniques should now be vindicated. Being 
able to write fairly explicit formulas for transfer functions and being 
able to see immediately the effects of feedback are compelling enough 
reasons to justify their use. 
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V. PROPERTIES OF THE ALGORITHM 


A . Stability 

Assumption 6 leads to stability in an easy way. Observing 
the resultant compensated plant leads to the observation that its 
characteristic equation's roots are: -a , roots of the model's character- 

istic equation, or contained in num det H(sl-F) . By the assumption, 
num det H(sl-F) ^G has all of its roots in the left half plane, hence 
the compensated plant is surely stable if the model is. The condition 
is not, however, necessary. For example, if p(s) equals a polynomial 
formed from all right half plane roots of num det H(sl-F) *G , and if 

i_i 

p(s) divides \ , we may cancel those terms in eq . (*) ,Fig.4 . 14 , 

i=l 

so that the roots will not appear in the result. 


B. 


Existence 


Assumption 5 guarantees that A^ exists and is non-singular. 


To show this note that: 


A i - H i Vi 


(5.1) 


-1 A H 

where det A(sI-$>) T = 'y 


s’*" , if H^ is non-singular. 


But : 


num det H(sl-F) 1 G 


= det H. 


s=0 


(5.2) 

\ 


num det A(sI-$) 1 r 


s=0 


= det A 1 = 


(5.3) 


Hence the assumption says that det H ^ 0 and det A^ ¥■ 0 . 
A^ exists and is non-singular. 


Thus 
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Again, the condition is not necessary. For example, if 
= 0 and X, = 0 we may still have success by cancelling the common 
factor of s in equation (*), Fig. 4.14, and renumbering (A becomes 

/ft 

A^ , A >2 becomes X and so forth). More complicated situations may 

also occur, but the essential ingredient is for T to be non-singular. 

£ 

This would have been quite difficult to state at the outset, hence we chose 
to use an intuitively appealing sufficient condition. 


Solution Invariance 


In an interesting case, the solution is largely invariant 
under changes in the plant. Recall Fig. 4.16, which shows the compensator 
as a transfer function from u and x to u . Suppose that the plant 
and model are both given in equicontrollable forms of the same state 
dimension. Then if we had written the compensator as: 


■w: 


u 


• -m: 


w 


and the plant and model as: 


PLANT: [ H ] 


MODEL: [ A ] , 

where q is a vector of parameters, then: 


Bf BG Bh 

5^7 = 9 ’ 5^7 = 9 ■ 557 


"0 t I “ 
i 
i 


0 


i 


. F(q) _ 


Li J 

”0 . I “ 
1 


' 0 

_ 

t 




_ i . 


Bj 

c _ c 

* ’*T 


Bf (q) 
~Sq~" 


(5.4) 


(5.5) 


(5.6) 


since A^ , J , etc., depend only on H and not on F , and T Q = I . 
Moreover, K was partitioned so that K 1 was nxm ! Hence K 2 vanishes. 
Therefore, the reduced compensator elements F ,G ,H are invariant under 
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changes in q , whereas J c is affected linearly. Note that we may 
relax the above restrictions to the point where: 

1. The plant is equicontrollable, and of the form: 



2. The original model's dimension is less than or equal 
to that of the plant. 

The result of this observation is that one need only estimate accurately 
' the parameters of H . The parameters of F do not as strongly effect 
the form of the compensator. Since varies linearly in q. , we 

know that small errors in estimating the parameters of F will result 
in small errors in the compensator. Moreover, if better estimates of 
F become available it is a trivial matter to adjust the compensator 
to match, since only the "feed- through" gains J c need be adjusted. 

D. Comparison with the Classical Approach 

As discussed in Chapter I, the classical control problem, 
which led to the "modern" version we have considered, was not readily 
soluble. In fact, our algorithm cannot solve it either, since too much 
structure was placed on the solution. We shall, however, consider a 
compromise problem which we can solve, and which has many of the features 
of the classical problem. In itself, it is of some interest since it 
includes a reference input. Such a feature is common to many control 
problems. The block diagram of the problem we shall consider appears 
in Fig. 5.1; 



Fig. 5.1 

A Classical Problem 
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R(s) represents a linear system which generates the reference signal 
r(t) in response to impulses at a(t) . (a(t) can be thought to 
establish the initial conditions placed on the system) . Our problem 
is to make e(t) behave in some arbitrary manner by the choice of 
C(s) . For example, the designer might like e(t) = £ -1 [Q(s)] for 
some transform Q(s) . Suppose that is indeed the case. For 
notational purposes, we might say that we wish the transfer function 
from a to e to be Q(s) . Hence we wish to find C(s) so that if 
the modified transfer function from e to y were V(s) , then 

Q(s) = (I+V(s)) -1 R(s) (5.8) 

Hence, if we could make V(s) = R(s)Q 1 (s) - I , then the problem would 
be solved. So we have the model: V(s) = R(s)Q *(s) - I and the 

plant: W(s) . This is then a problem which falls under the ability 

of our algorithm. We need only check that the various assumptions and 
conditions of that algorithm are either satisfied or easily accommodated. 

A short study of the matrix Q(s) will be of some help. 

Q(s) represents the dynamic response which we would like the error to 
have. As such it is not unreasonable to suppose that an acceptable 
Q(s) might be diagonal (noting that, by the definition used above, 

Q(s) is indeed square). Moreover, a designer might well be pleased 
with a response of the form Q(s) = I/(s+a) . Although one usually 
would think of such a response as a bound on performance, we mean it 
here to be the actual dynamics of the error. 

Another point worth a brief discussion is the transform 
R(s) . The reference system is simply a model of the specified refer- 
ence input. In that case it is not unreasonable that R(s) also be 
diagonal: that is, that the separate input (reference) signals r^(t) 

could be separately described by the diagonal entries of R(s) . 

Although a(t) really only represents the initial conditions on the 
systems described by R(s) , it is not any problem to allow a(t) to 
represent any input. 

Putting these observations together, we find that: 
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or 


(5.9) 


V(s) = (s+a)R(s) - I 
V i (s) = (s+a)R (s) - 1 

If the state dimension of R^s) is greater than or equal to one, 
there is no problem of having a non-realizable model. This does 
not appear to be a problem, since it makes little sense for the 
reference input to be non-dynamic! 

In conclusion then, it is quite possible to wrestle the 
classical type of reference input control problem into a form which we 
may easily handle. The algorithm is used to solve an auxiliary 
problem which, in turn, makes the main problem trivial. Note that the 
nature of the problem involves a decoupling (since the defined model 
is diagonal). This points out the subserviance of the decoupling 
problem to that of model following, (see Appendix C) . 

When there is no reference input, the problem is amenable 
to our algorithm. The classical approach seems, at first, to be simpler 
if for no other reason than the compensator is easy to define. On the 
other hand, the indicated inverses are especially tedious and the 
resultant system is usually of over-large state dimension. 

As an example of how the result of using the classical 
approach compares with our design, consider the following problem. 




- 


1 

s+4 

W(s) = 

s + 

5 (s+1) (s+2.5) 


n 

s+4 


u 

(s+1) (s+2.5) 


1 

2 

V(s) = 

s + 

1 s + 2 


0 

s+4 


(s+2) (s+3) 

Using the procedure shown in 

Chapter I, we would find 


\ 
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C(s) = W(s) _1 V(s)[ I-W(s) ] 

s + 5 
s 

0 

This solution requires 3 states. 

On the other hand, using our algorithm, we would find that 
only one dynamic element, namely an observer root, is required. If 
0(s) is an observer, then our solution is shown in Fig. 5.2; 


( s+5)(s +4s+6 ) 
s ( s 2 +4s+2) 

( s+l)(s+2.5 ) 

(s 2 +4s+2) 


(5.11) 



Fig. 5.2 

Form of Solution to the 
Classical Problem 


where C is defined by: 


1 

1 1 4 -1/3 

i 

1/3 

A “ 

U 

0 

i 

1 ! 0 -1/3 

1/6 _ 

/V 

X _ 


(5.12) 


This feedback law was generated by our program using suitable realiza- 
tions of V(s) and W(s) . 

It is worth noting that our solution is much simpler than 
the classical design. A quadratic loss approach would result in a 4 
state compensator (3 for the model realization, and 1 for an observer). 
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VI. COMPUTER PROGRAM AND EXAMPLES 


A. The Computer Program ^ 

The foregoing algorithm has been programmed for an IBM 
360/67 at Stanford University, in Fortran IV. The details of the 
program are hardly germain here - particularly in view of the brute 
force nature of the programming. Besides the computation of the con- 
troller, the program includes an impulse response simulation for the 
model and compensated plant. This provides a graphic check as to how 
well the solution actually matches the model's behavior. Some controls 
on roundoff error are included such as using double precision throughout, 
setting "canonical elements" in a Luenberger form to 1 or 0 as the case 
may be, and occasionally rounding off matrices to set small elements 
to zero. 

A flow diagram is not particularly useful in describing 
the program since only the reduction to Luenberger form, a subroutine 
to find the independent columns of a matrix, and the simulation routine 
are iterative in any sense. Hence the flow is virtually always down such 
a chart. The program was a literal translation from the description of 
Chapter IV to Fortran. Subroutines were used heavily. The following is 
a rough scenario which is a reverse translation from the main program 
back to verbage. 

1. Perform the requisite input, immediately echoing back 
most of what was read in as a check. 

2. Reduce the plant and model to Luenberger canonical 
form, [13], (It turns out to be easier to go 
through Luenberger form rather than to equicontrollable 
form directly - see Theorem 1, Chapter III.) 
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3 . 


Print the results. 

Add states as required and get the plant and model 
into equicontrollable forms of the same state dimension. 

4. Perform a roundoff to set small entries to zero. 

Print the result. 

5. Compute the polynomial matrix: 


ladj 




- ■ ' * 

£ „ i-1 

\ 

£ i_i 

\ H. s 

I* 

\ A. s 

Z i 

I 

Z 1 

i=l 

/ 

i=! 


A V- a s i-l 


2 v 


i=l 


and the corresponding 


det 


2 
i =i 




(Note: k is determined and adjusted as per Chapter 

IV, section C) . 

6. Again perform a roundoff. 


7. Check for factors of s and stability of 


P 

2 

i=l 


X. s 

x 


i-1 


If assumptions fail — exit the program with error message. 

8. Print A. and X. ■ . 

9. Reduce model to bar-star coordinates. 

10. Invert A^ . If it is singular — exit with message. 

11. Compute compensator. (This is a vast understatement. 

The code is lengthy, the essence is bookkeeping. 

Print the result. 

12. Minimize the dimension of the compensator. Print it. 

13. Find the compensated plant and simulate it. 

Further details of the program may be gleaned from consulting 
the Appendix B. It gives a description of many of the matrices and 
variables computed. 
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B. 

1. 

problem. 


Examples 
Example 1. 

Let us use the algorithm described above on a fairly simple 
PLANT: 

F 
H 



MODEL: 


$ = 


A = 




The first step is to reduce these to a Luenberger form. 

1 


PLANT 


f l 


h l = 


MODEL : 


-5 

0 


0 0 
0 1 
0 2.5 -3.5. 

14 1 

0 4 1 


*L = 


^L 


-1 

0 

. 0 
1 
0 


0 

0 


0 

1 


-6 -5 J 

6 2 
4 1 


g l = 


1 

0 

0 


r L = 


1 

0 

L0 


O' 

0 

1 


0 

0 

1 J 


V 

( 6 . 2 ) 


Now we add a root at s = -15 to each system (a = -15) and again 
reduce to a canonical form similar to Luenberger * s . 
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PLANT 2 : 


F = 


-20 0 0 -75 

0 0 10 

0 -2.5 -3.5 0 


10 0 0 



0 

0 

1 

0 




4 

4 


1 15 

1 0 


-16 0 0 -15 




(6.3) 


A = 



6 

4 


2 15 

1 0 


A simple state renumbering will reduce the above to equicontrollable 
form. Note that only one state was required since the previous block 
lengths of 1 and 2 (see Theorem 1, Chapter 3) indicate that 4 states 
will be needed. The equicontrollable form is then: 


PLANT 3 : 


ro 



0 1 

0 0 

0 -20 


0 

1 

0 


0 -2.5 0 -3.5 



0 

0 

0 

1 


H 


* 


15 

0 


4 

4 


1 1 
0 1 


(6.4) 


MODEL 


3 



> 


0 0 
0 0 
1 0 
0 1 


A 


* 


15 

0 


6 

4 


1 2 
0 1 
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Up to this point we have simply been maneuvering the plant and model 
into forms which are convenient, while maintaining their transfer 
functions. Now we proceed to find the compensation required. It is 
easy to compute that: 


A 3 = 

1 

0 

1 

f 

1 

> 

to 

li 

19 

0 

6 

19 

\ = 

o o 

CD 

1 1 

8 ' 
60 

X 3 " 

1 

> 

li 

19 

} 


60 



is clearly invertible, hence there will be no difficulty with 
getting the compensator. We can see that there are no zeros in the 
RHP in num det W(s) = s + 3 , so stability of the solution is sure. 
Now via Fig. 16 we may simply write down the solution, after checking 
the definition of some assorted matrices. To wit: 


J 


11 


[ 0 ] 


J 12 = [1 0 0 ] 


( 6 . 6 ) 



1 

0 ! 

0 

0 


’ 0 

o ; 

i i 

0 

i ; 

0 

0_ 

J 22 

0 

0 

i 

0 . 
o ! 

0 

0 1 






0 

o ; 

0 1 


A simple observation aids in writing the solution. Consider: 



Ik 

+ a i 

{ t" 1 ] 

, 1 ! (A, K 

+ Jo, 1 


A 2 

k|_ 

1 1 J- 

2 / : V k 2 

21 ) 


By computation (hand or machine) we find that H c = [0] . Thus only 
feedback will be required (along with a suitable input transformation). 
The solution is then: 


~1 

1 ; 4 -2/3 

i 

2/3 " 

_0 

1 | 0 -2/3 

1/6. 


( 6 . 8 ) 


written in original coordinates. One may check that indeed this does 
the job. The above matrices were all taken from our computer program 
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compiled by a load-and-go type processor, and run on an IBM 360/67. 
The execution took 1.34 seconds. (Compilation took 8.00 seconds). 
Execution time would have been of the order of 0.3 seconds if it had 
been compiled using the IBM Fortran IV Level H compiler-link editor, 
but compilation would have taken over a minute. 


2 . 


Example 2 : 

This example first appeared in a paper by Tuler and Tuteur 


[19] 

. The problem is 

defined 

by the 

equations: 



^ 0.0 

1.0 

0.0 

0.0 



0.0 

-2.93 

-4.75 

0.78 


F = 

0.086 

-0.0 

-0.11 

-1.0 



_ 0.0 

-0.040 

2.59 

-0.39 



0.0 

0.0 





0.0 

-3.91 




G = 

0.035 

0.0 





-2.53 

0.31. 




H = 

‘ 0.0 

1.0 

0.0 

0.0 



_ 0.0 

0.0 

0.0 

1.0 



0.0 

1.0 

0.0 

0.0 


$ = 

0.0 

-1.0 

-73.14 

3.18 



0.086 

0.0 

-0.11 

-1.0 



0.0086 

0.086 

8.95 

-0.49 



’ 0.0 

0.0 





0.0 

-3.91 




r = 

0.035 

0.0 





-2.53 

0.31 





■ 0.0 

1.0 

0.0 

0.0 


A = 

. 0.0 

0.0 

0.0 

1.0 

The 

solution we obtain by use 

of our 

■ program 

is: 


z = 


u = 


F z + G 
c c 


X 

/S 

u 


n * + o 
c c 


(6.9) 


( 6 . 10 ) 
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I 


where: 


F 

c 

G 

c 


H 

c 


[- 

[ 


24. 92975 J 
3.8066 034 


-0.4247245 j 3.142755 
-18.09820 


-2.441089 

36.66145 


-0.2604497 
-36.66145 J 

( 6 . 11 ) 


1.000000 

0.000000 


0.00000000 J 0.3074326 

i 

1.00000000 1-4.6682222 
-0.9029511 
25.48522 


-0.1369023 

-0.09383566 

-3.715906 

54.65753 


If F ,G ,H J are used to double precision (as they are computed, but 
not as they are listed above) then the model and the compensated plant 
agree to 14 decimal places in their respective impulse responses (over 
the first 4 time constants of the transient). Surely this is a good 
fit. Moreover the solution was obtained after one pass of the program 
(nee; algorithm), and took but 6 seconds of 360 time. The solution 
put forward by Tyler and Tuteur was only obtained after some unspecified 
number of man-machine iterations (i.e. new choices of the Q and R 
matrices), and involved 4 states in the compensator (by virtue of the 
quadratic loss approach). The observant reader might notice that no 
output matrices were named there. On the other hand, the solution they 
obtain, which purports to be matching all of the states, actually 
succeeds mainly in closely matching the characteristic equations rather 
than the transient responses, much less transfer functions. Our 
solution above, matches the first three states nearly perfectly in a 
transfer function sense. Should it be required to match the third 
state, say, instead of the fourth, one would simply make; 


0.0 0.0 0.0 0.0 

0.0 0.0 1.0 0.0 


( 6 . 12 ) 


This may then be computed to have the solution: 
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F 

c 


G 

c 



|^-72.62536 j 

£ 0.03243403 

0.08903046 

’ -2.771866 

-0.07890249 

1.00000000 

0.00000000 

-0.1098584 

-0.4967333 


0.150421 

2.771866 


0.00000000 
1 . 00000000 
-0.3698480 
17.48052 


-0.0001956313 
-0.1186034 1 


-0.00295 

-0.00008405342 

-0.04065118 

0.6149679 


(6.13) 


This solution is not surprisingly quite different from the first. 
However, one may as easily solve for any compromise between the two 
or any other desired output set - so long as it is limited to two 
directly. (Since one state is the integral of the other, it will be 
matched up automatically if its associate is). 


3. 


Example 3 ; 

Rynaski and Whitbeck [17] present a similar problem. 


F = 


G = 


H = 


$ = 


-0.751 

0.0000046 

0.000572 

-1.604 

1.0 

0.0 


0.0 

0.0 

0.0 

32.2 


-0.0296 

17.45 

1.0 

-0.0000214 

-0.0009599 

-0.681 


- 




0.0015 

-2.65 




0.0 

0.0 




8.46 

0.0 




-0.0069 

-0.0326 _ 




r i.o 

0.0 


0.0 

0.0 "1 

L o.o 

0.0 


1.0 

0.0 J 

-1.28 

0.0000044 

0.0005067 

-0.2558 ’ 

1.0 

0.0 


0.0 

0.0 

0.0 

32.2 


-0.0263 

30.58 

1.0 

-0.0000086 

-0.0009859 

-0.8584 


To wit: 


(6.14) 


/Contd . 
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/ Contd. ..(6.14) 



A = 


[ 


0.0048 

-0.898 



0.0 

0.0 



17.4 

0.0 



-0.00934 

-0.129 



1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 

0.0 


Again, the answer yields remarkably good results, namely 7 decimal 
places accuracy as measured by the simulation over a similar time 
period as the previous example. The output matrices had to be contrived, 
but exactly the same arguments may be put forward as for Example 2. 

Also it is possible to match 3 of the 4 states nearly exactly, and there 
is considerable freedom as to which states or combination of states will 
be matched. The solution obtained is: 


-0.8044424 ^J 



[-0.0115643 

0.0751778 

| -0.8354341 

0.0443522 

0.0008648 

0.5827508 

-0.8354341*] 



-0. 0227829 J 



2.0587376 

0.0000000 

j 0.4112750 

-0.0006471 

0.3388679 

j 0.2127252 

0.0424707 

0.0002993 

1.0295431 

0.0011583 

0.0000224 

-0.5221243 


(6.15) 



APPENDIX A 

A MINIMUM NORM APPROACH 


A. The linear problem . 

Given the plant and model: 

PLANT : x = Fx + Gu 

MODEL: £ = <J>£ + G . 

We might wish to find a feedback control law u = Kx 3 : 

T II OT (F+GIDt.,2 1|aN 2 r .T . 

J = ||e - e ||q ; ||A|| = trace [A QA] 


(6.16) 


(6.17) 


is minimized over choices of K . This is a generalization of Erzberger's 
[ 6 ] approach. We refer the reader to his paper for motivation of the 
formulation, it incorporates two distinctive features not possessed by 
hi s , however , 

1. A weighting matrix Q is included (Q ]> 0) so 
that one's relative interest in matching the states 
can be expressed. 

2. A parameter t is included which is a gross 
measure of closeness. As t -» 0 , this problem becomes 
the same as Erzberger's. 

The main hurdle to be overcome in solving this problem 
is finding the gradient of the cost function with respect to K . 

Further we would like the gradient to be linear in K . To that end, 
let us make a model of the plant. 


x(n + 1) = e FT x(n) + F 1 (e FT - I)Gu(n) (6.18) 

(Note: F -1 need not exist for the term to exist, as a series 

expansion will confirm). 
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Suppose we choose t to be small so that the above 
is a good representation and further that the feedback law: 


u(n) = Kx(n) + u(n) 


(6.19) 


is desired. Then: 


x(n + 1) 


Ft -1 Ft _i 

= e + F (e - I )GK x(n) + F 


+ F _1 (e FT - I)Gu(n) 


( 6 . 20 ) 


u(n) represents external inputs. Now, the continuous system, after 
such feedback, is: 

x = (F + GK)x + Gu (6.21) 

But this system is virtually equivalent to the preceding (for small t), 
hence: 


FT -1, FT [F+GKlT 

e + F (e - I)G2-e L J 


( 6 . 22 ) 


Then: 


J= He $T - e FT - F“ 1 (e FT - 1)GK||| = ||e* T - e*^ 
-2<r T Q(e <5t - e Ft ) , K > + <r T QTK, K> 


(6.23) 


■1, Ft 


where f = F "(e* I)GK and <(A,B^ = Vtrace A T QB 


Then taking V^(J) formally , we obtain: 
K 


V Tr (J) = -2r T Q(e 0T - E FT ) + 2r T QTK 


K 


Setting Aj^(J) =9 we obtain: 


K = [rV] e F ] 


(6.24) 


(6.25) 


This is then a solution to a form of model following problem. 
It suffers from the basic difficulty of Erzberger's approach, namely 
that the plant and model must have vast structural similarity for it to 
be successful. Fortunately, however, there are many problems having 
this trait. In aircraft control problems we have a great deal of 
information as to structure. One often finds the plant and model are 
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different only in a few parameters, while they have the same essential 
organization of dynamics. Hence the success of Tyler's problem (among 
others) is not surprising. 

One might notice that, as in the quadratic loss formulation, 
we have an unknown weighting matrix Q . No attempt will be made to 
minimize the problems engendered by it. Suffice to say that the solution 
here takes less than one second to perform, hence a time-share mode 
might be quite feasible, wherein the designer would rapidly cut and try 
various Q's until an acceptable solution appeared. That such an 
approach can yield spectacular results will be shown by example. 

Rynaski, et al [16] give an interesting problem, depicted in 
Pig. A.l. They give the usual solution, but indicate that improvement 
in accuracy would involve substantial increase in control effort. 

However, trying our approach yields an interesting result. Not only 
can we find a pure feedback control which does slightly better than 
theirs, but it does so with uniformly less control effort. For an 
increase of about 20$ in maximum effort we may do a spectacularly close 
match. Clearly we may choose any tradeoff in between that we desire. 
These results are depicted in Figs. A. 2 and A. 3. Plotted are the pitch 
responses of the open loop aircraft, the model, Rynaski 's solution, and 
the two solutions obtained via our method. The associated control 
efforts are plotted to verify our claim. 
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Fig. A.l 
Rynaski Problem 
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A final example, from Rynaski and Whitbeck [17], is almost 


trivial . 


F = 


$ = 


-0.751 

0.0 

0.0 

-1.60 

1.0 

0.0 

0.0 

0.0 

0.0 

-32.2 

-0.029 

17.45 

1.0 

0.0 

0.0 

1 

O 

0) 

00 

-1.285 

0.0 

0.0 

-0.26 

1.0 

0.0 

0.0 

0.0 

0.0 

-32.2 

-0.026 

30.6 

1.0 

0.0 

0.0 

-0.86 


G = 


with Q = 


0.0015 

0 

8.5 

-0.0069 

1 0 
0 0 
0 0 
0 0 


was 


-2.65 

0 

0 

-0.033 

0 0 
0 0 
1 0 
0 0 


(6.26) 


the resultant system 



-1.275 

0.0 

0.0 

-0.264 

GK = 

1.0 

0.0 

0.0 

0.0 

0.066 

-32.2 

-0.0260 

30.53 


0.993 

0.0 

0.0 

-0 . 674_ 


(6.27) 


Enough said.' 


In conclusion then, if a problem with such a high degree of 
structural consistence between plant and model appears, this method is 
worth trying. The computation times involved are so short that it might 
be used as first cut in any event. If five or so tries fail to yield 
an acceptable solution, it is perhaps best to switch to the more general 
algorithm. 


93 - 



B. 


A non-linear extension 


Suppose we have the plant: 
x = Ff (x ) + Gu 


where where 


w 


f(x) = f 2 (x 2 ) 


f (x ) 
n n 


(6.28) 


and F and G are the usual matrices, and a model: 
where 


S = ®f(S) 


f i ( ^i ) 


Jet) = 


f (S ) 
m m 


m ^ n 


and further assume that for the output matrix H , 
f (Hx) = Hf (x) 

Then say we desire £(t) = Hx(t) given that |(0) = Hx(0) 
implies: 

S(t) = Hx(t) 


But : 


= HFf (x) + HG u 

S = ®f(S> 

= $£(Hx) 

= <&Hf(x) 


So that: 

HFf (x) + HGu = 5>Hf(x) 


or 


define 


HGu = (OH - HF)f(x) 
u = Kf(Dx) 


(6.29) 


(6.30) 


(6.31) 


(6.32) 

(6.33) 
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where 


/s 

A 

f(Dx) = 


f ^Ci-jX) 


f (d x) 
P P 


= Df(x) 


(we assume such a D exists). 

Then we desire a K 3: 

[HGKD - <4>H - HF)]f(x) = 0 for any x . 
It is sufficient that: 

HGKD = 3>H - HF 

Then the solution, if it exists is: 

K = (HG) + (4>H - HF)D + 

For example, the systems: 


- 

'-1 1 - 

x. 

— 

0 

x. 



l 


1 

= ■ 


2 

+ 

u f(x) = _ 


M 

O 
l 

X 2 J 

1 

X 2 J 



"-1 

1 " 

h " 

= 



2 


3 

- 

0 

^2 J 


H = 


1 

0 


so f(Hx) = f(x) 


and D = I is also reasonable. So K = 


2 

[2 


0 

0 . 
0 ] 


So u = [2 0] f (Dx) = [2 0] 


2 

X 2 


2x, 


(6.34) 


(6.35) 


(6.36) 


(6.37) 


(6.38) 


(6.39) 


(6.40) 
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which is correct, although trivial. Nonetheless, the idea is that in 
cases in which the previous linear approach is used, we may allow some 
non-linearity into the problem without changing the solution. Such a 
condition is often present in aircraft problems; i.e. a small amount 
of nonlinearity. This might be effectively modeled in the form we 
suggest which would assure us that if the linear problem could be 
solved, then the non-linear one could be also. 



APPENDIX B 


DEFINITION OF SOME VARIABLES 


Our purpose here is to describe those variables which are 
used in the algorithm which are particularly obtuse in definition. 
They will be considered in order of occurrence in the original. 


( v 


The number of states added in z depends on: 

1. The number of states n^ needed to make m|n+n^. 


n^ is chosen. 


The smallest such 
2. Whether or not the model has more than n+n 


n more, then n 
2 2 


extra states 


states. If it does, say 
are added. 

We suppose that the system is in Luenberger form and that a 
state map is available. It has p rows and m columns 
(p = (n+n +n )/m). It represents how the current block 
diagram of the plant is numbered. For example: 


4 

V 

zL 

i=l 


a .x . +u„ 
i i ] 


*± 

1 

i=l 


b.x. +u 

li 2 



This would yield the state map: 

‘ 1 4 “ 

2 0 

.3 0 . 

The zeros are merely place keepers. To illustrate how states 
are added, and in turn how is created, suppose in the 

above example we wish to add 4 states. The new map would be: 
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5 


6 


The procedure is to "push down" the short blocks and to add 
states into the spaces required, numbering left to right and 
top to bottom. The block diagram then becomes: 


Ivv 

i=l 

2 vr 


1 

X 3 

1 

X 2 

1 "l 

1 

s 


s 


s 

s 


1 

X 4 

1 

X 8 

1 

X 7 1 

s 


s 

— • — ■* 

s 

— s 


This is reminiscent of the discussion after Theorem 1, 
Chapter III, we may define: 


■ J 12! J n] 


The partition is made between the n^ and (n+l) st 


columns of J, 


Even though no a has been mentioned, 


the dynamics associated are contained in a^ and b^ rather 
than in . In actuality the added states would have been 
l/Cs+a) rather than 1/s , but it is an easy matter to separate 
those dynamics from the structural information of (which 

is invariant under a) . 

is entirely analogous for the model. The notational 
changes should be apparent. 
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^ >1 


e i e . . . e 

o e 1 - -- 0 


F F 
1 2 


then 


r = - f 

i L i 


F 9 ”• F 

2 p 


KT^ is defined rather than K for notational convenience. 



This is subject to the restrictions: 

1. A,, = 1 

k 

2. r is to be as small as possible. 

3. The left side of the above is to be minimal. That 
is, any cancellations possible are made between 
numerator and denominator. 

4. The numerator and denominator on the left are to be 
the same length (as shown) . 
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X.I 

1 



k-1 blocks 


all blocks are mxm 


k-1 blocks square 
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p + k - 1 blocks square 

Note that T is nonsingular iff A is nonsingular. Also 

Ci 

T is (p+k~l)m dimensional . 
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APPENDIX C 


DECOUPLING 

After becoming familiar with the technique of this 
dissertation, it becomes clear that the decoupling problem (see 
Wonham and Morse [22], and Morse and Wonham [15]) is closely allied 
to model-following. To wit, the model is decoupled. Granted the 
decoupling problem specifies no particular model. On the other hand, 
we can easily see that if the plant satisfies the conditions given in 
Chapter II, then virtually any decoupled model can be followed perfectly. 
This observation is borne out in theory. In [22] a necessary and 
sufficient condition is given such that a system may be decoupled by 
dynamic compensation. We have shown (see [23]) that their condition 
is equivalent to ours under the input and output restrictions given in 
Chapter II. Hence our intuition is nicely fulfilled in reality. 

Specifically we may state the following theorem, 

proven in [23]: 

Theorem : Given the system (H,F,G) with H mxn , 

F nxn , G nxm and full rank, then (H,F,G) may be decoupled by a 
dynamic compensator C(s) , of the form given in Fig. 2.1, if and only 
if num det H(sl - F) 1 G has none of its roots at the origin of the 
s-plane. 

This theorem constitutes a simple test for decoupling 
in the m-input, m-output case. The design is as per the previous 
solution (Chapter IV) with the model virtually arbitrary. Note further 
that the comments in Chapter V also may be applied to decoupling. In 
particular, Wonham and Morse and the above theorem do not specify 
that the resultant decoupling be stable. If that is required, the 
further restriction of the roots of num det H(sl - F) ^G to the left 
half plane will suffice. 
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